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An analytic means of studying the resonant decay of the inflaton field is developed for the 
case of background expansion, H 0. It is shown that the parametric resonance in the 
inflaton’s decay need not disappear when the expansion of the universe is taken into account, 
although the total number of particles produced is fewer than in the H ~ 0 case. 
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1 Introduction 


Soon after the first models of inflation were published in 1981-1982 0 §, several papers were 
written to study how the universe would reheat following the supercooled inflationary state. [|3[] 
The process appeared to be straightforward: in models like new inflation (and chaotic inflation Q) 
which incorporate a second order phase transition to end inflation, the inflaton field would wind 
up oscillating around the minimum of its potential near the end of inflation. These oscillations 
would produce a sea of relativistic particles, if one added (by hand) interaction terms between the 
inflaton and these lighter species. The decay rates could then be calculated, where <fi is the 

(decaying) inflaton field, x is some light boson field, and is some light fermion field. From these 
decay rates, the energy density ( p p ) of the produced particles could be calculated, and related to 
the final reheat temperature (T r h): p p — vr 2 A ? e ff T^/30, where !V e fj- counts the effective number of 
massless spin degrees of freedom. This method was thought to estimate T r h to within an order of 
magnitude. 

Recently, much attention has been focused on an overlooked but dramatic feature of the 
production-via-oscillation model: the inflaton’s oscillations should be unstable, and should exhibit 
parametric resonance. Certain modes Xk within small bands e <C k should grow exponentially, 
swamping the production in other modes. This effect has been studied by Kofman, Linde, and 
Starobinsky in || for both the broad and narrow resonance regimes, and, following different meth¬ 
ods, by Shtanov, Traschen, and Brandenberger [p] for the narrow resonance regime. The importance 
of “induced amplification” in the inflaton’s decay has been addressed recently with still different 
methods by Boyanovsky et al. in pi] [1] . See also [9] flOl] . Earlier work on some of these questions 


may be found in [11] [|l^]. The crucial feature of the new reheating scenario is that the equation of 
motion for the light boson field will not obey the simple harmonic oscillator equation. Instead, for 


an interaction of the form C' = —\g 2 (t> 2 X 1 , the mode functions of the boson field will obey: 


Xk + 3-ffxfc + 


' k 2 


+ g 2 4> 2 Xk = o, 


(i) 
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where overdots denote derivatives with respect to cosmic time t, a(t) is the scale factor for the 
(assumed) flat Friedmann-Robertson-Walker spacetime, and H = a/a is the Hubble parameter. 
Since the inflaton field ( 4>(t )) is rapidly oscillating, for frequencies loj H the 3 H\k term may be 
neglected and a(t ) treated as constant. Then this equation for the boson field mode functions will 
look like the well-known Mathieu equation (see, e.g., [|b|] fl4}| ): 

Xk + (1 + h cos cojt) Xk = 0. (2) 


As demonstrated in i i jIT| and treated below, for loj = 2u>h + e, with e <C u>k, solutions grow 
like Xk(t ) oc exp(±[i±t), instead of like Xk(t) oc cos(u>kt). This means that the number of x bosons 
produced goes like N™ s oc exp(2/i + t), so that the decay rate for <fi —> 2y is much larger than 
originally calculated: T^ es T x . (Note that this parametric resonance is only effective for inflaton 
decay into bosons; the exponential increase in fermion modes is forbidden because of Fermi-Dirac 
statistics.) This large decay rate increases p p , ultimately giving T r r ^ s 2> T r h- 

In this paper, we develop an analytic means of studying the inflaton’s decay into other particles 
for the case of H ^ 0. We find that the total number of particles produced from the parametric 
resonance effect when H ^ 0 is fewer than in the H ~ 0 case, but can still be exponentially greater 
than when the resonance is neglected altogether. Rather than expansion of the universe, then, the 
true threat to resonant inflaton decays appears to be back-reaction of the produced particles on the 
amplitude of the decaying inflaton field. While we do not study this back-reaction effect in detail 
here, some approximate schemes for including some of its effects are discussed in the Appendix. 

We concentrate on the case of the inflaton field cj) decaying into further inflaton bosons, via 
the self-interaction potential V(4>). We restrict attention to the small-amplitude, narrow-resonance 
regime, as the inflaton field oscillates near the minimum of its potential, similar to that studied 
in HJ. In section 2, the number of particles produced per mode k, Nf., is calculated for non-zero 
expansion, assuming that the decay is not resonant. Section 3 presents the calculation of N^ es 
when the expansion of the universe is neglected. In section 4 we calculate N^ es taking into account 
the expansion of the universe, and demonstrate that the parametric resonance need not vanish. 
These values of correspond to particle production during the narrow-resonance regime of the 
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“preheating” epoch discussed in [j5] [15]. Final reheating temperatures can then be calculated, as 
discussed in ||, using the methods developed in [i|], where these methods should now be applied 
to the produced particles of the preheating epoch; this will not be included here. For more on 
the question of thermalization and the final reheat temperature, see also ||. Concluding remarks 
follow in section 5. 


2 Non-Resonant Inflaton Decay in an Expanding Background 


We will study a very simple model of inflation: chaotic inflation with the potential V((j)) = |A</> 4 . 
We will assume that the metric may be written in the form of a flat Friedmann-Robertson-Walker 
line element, ds 2 = g^dx^dx 1 ' = —dt 2 + a 2 (t)dx 2 . Then, given the lagrangian density 


£ = 


1 

167tG 



v(d>) 


(3) 


the equations of motion take the familiar form: 


H 2 

4 > 


87 tG 


V{4 o + 


1 


dV 


+ 3 H(f> - <j) H—j— — 0. 

or dcp 


(4) 


Next we decompose the inflaton field into the sum of a classical background field and a quan¬ 
tum fluctuation: 4>(x,t) = ip(t) + 8cj)(x,t). (We will neglect the metric fluctuations which couple 
to these inflaton fluctuations for the present analysis; for more on these metric fluctuations, see, 
e.g., [|l(J |l7]].) In order to study the production of <fi particles when the inflaton field is oscillating 
near the minimum of its potential, we will need to calculate mode solutions for the quantum fluctu¬ 
ations 5(f> before and after the phase transition, matching the solutions by means of a Bogolyubov 
transformation. 

We will first concentrate on the modes for the fluctuations near the end of the inflationary 
period. In terms of conformal time, dg = a~ 1 dt 1 the fluctuations are quantized as follows: 


5^>(x) 


d 3 k 

(27 t ) 3 / 2 


S(j>k(V) H e 


^ + S</>Uv) a\e~ ik -* 


(5) 
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where hats denote quantum operators, and the creation and annihilation operators obey the canon¬ 
ical commutation relations: 


a k > a i 


at, at 
k ’ i\ 


= 0 , 


dj^ , d 


= 5 3 (k — £) , with ar|0 >= 0. 


( 6 ) 


From equation (||), and denoting d/dr] by a prime, the equation of motion for the mode functions 
54>k{v) becomes: 

8(f)I + 2— Sep/ + f k 2 + 3A(/3 2 a 2S ) 5<pk = 0. (7) 

a V / 

Next we introduce a “conformal field” ip = a5<p, whose mode functions obey the equation of motion: 


V’fc + 


k 2 — ——h 3A tp 2 a 2 
a 


ipk = o. 


( 8 ) 


Near the end of inflation, the scale factor a[t ) will not in general retain its de Sitter (exponential) 
form. Instead we may write a(t) oc t p , with 1/2 < p < oo, where p = 1/2 corresponds to the 
radiation-dominated epoch, and p —> oo recovers the de Sitter epoch. This scale factor corresponds 
to a conformal-time scale factor of a{rf) oc rf Because we are ignoring the parametric reso¬ 

nance from the classical (^-field’s oscillations in this section, we may simply assume that near the 
minimum of its potential, p 2 <C k 2 . If we lastly define a new field x = p~ l ^ 2 ipi and define a new 
variable z = krj, then equation @ becomes 

(3 p - l) 2 


d 2 Xk 1 <bCk 

dz 2 z dz \ z 2 4 (p — l) 2 


Xk - o. 


(9) 


This is now in the form of Bessel’s equation. Mode functions for the original field 5<p may then be 
written in terms of Hankel functions: 


Hk{ri) = 


V 


!/2 r 


V = 


a (v) l 

3p — 1 


A k Hi l \k V ) + B k Hl 2 \k V ) 


( 10 ) 


2(p — 1)' 

We assume that during inflation, the inflaton field is in its vacuum state (see section 5.2 of ]ll 
for further discussion of this point). The coefficients A k and B k thus may be fixed by choosing an 
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appropriate quantum vacuum state. The Bunch-Davies, or “adiabatic”, vacuum requires that the 
fluctuations ipkiv) behave as Minkowski-spacetime mode functions far inside the horizon, that is, 
ipk —> (2 &) 1 / 2 exp(— ikrj) for kp 1. [17] 18] [20] From the asymptotic properties of the Hankel 
functions, this gives 


7T 


Ak = 0 , Bk = exp 


7T 

—i— I v + 


1 


( 11 ) 


For the problem of (non-resonant) reheating, we will be interested in the opposite asymptotic limit: 
hr] -C 1, corresponding to long-wavelength modes. In this limit, the fluctuations 'ipkiv) take the 
form 


- 0 k(v ) 


2 v ~ l 

r ■ 7T 

( i\i 

exp 

VvF 

1 2 



r(i/) 77 1 / 2 {kifj v . 


( 12 ) 


In order to calculate the number (Nk) of (/>-particles produced per mode k, we match this solution 
to long-wavelength mode solutions following the phase transition into the radiation-dominated era. 

To do this, we make use of a Bogolyubov transformation. (See, e.g., |l^].) For long-wavelength 
modes, we may approximate the phase transition as instantaneous, occuring at some time 77 * [ 21 ]: 

77 < ? 7 * : a(rj) = (a 0 ri) p ^ 1 ^ 

V>V* : a(rj) = C (77 - 77 ), (13) 


where fj = 77 * — (a 2 ??*) -1 . It is convenient to set a(? 7 *) = 1, which sets 77 * = a“, and thus 77 = 0 
and C = a Q . (Note that in the limit p —► 00 , the reference scale a Q becomes —H 0 , where H a is 
the Hubble constant of the de Sitter spacetime.) Next we perform a Bogolyubov transformation to 
match the long-wavelength mode solutions in the radiation-dominated era: 


'iPk(r) > 77 *) = 


1 


V2k L 


a k e 


—ikrj 




-\-ikij 


(14) 


The Bogolyubov coefficients a*, and f3k may be determined by requiring that both ipk and i/)' k be 
continuous at 77 = 77 *. Then, using Nk = |/3fc| 2 , this gives 

4 V_3 / 2 „ 0 , . / 1 \ 2 fk ^~ 2u ~ 1 


N t = IAI 2 = 24—r 2 (0 (*- 1) (4) 


(15) 


Equation (15) is the main result of this section, and can now be compared with the cases in which 
the effects of the inflaton’s parametric resonance are included. 
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3 Particle creation from Parametric Resonance, with H ~ 0 


We return to equation (|j) for the fluctuations ipk(r]). Rather than neglect the 3\ip 2 (rfci 2 (rf term 
in this section, however, we study the effects of this term as the inflaton field ip oscillates near the 
minimum of its potential. If the frequency of these oscillations is large enough, then we may expand 
the fields ip(rf and if 2 (if) in conformal-time harmonics, analogous to the cosmic-time harmonic 
decomposition adopted in fh]. Keeping the lowest term, we may write 

ip 2 (rf ~ ip 2 cos( 7 ? 7 ), (16) 


where ip 2 is a slowly-decreasing amplitude. This quasi-periodic approximation for ip 2 (if may be 
compared with the numerical results (for Minkowski-spacetime) calculated in J?]]. For this section, 
we will assume that ip 2 — constant. (As emphasized in j|] [Q], the chief contributor to the decrease 
of ip 2 over time is the back-reaction from created particles, and will be further discussed in the 
Appendix.) Substituting this ansatz into equation (pi) yields: 


■ 0 k + wl 1 + ga 2 cos ( 777 ) - 


f>k = 0 , 


where 


2 _ j 2 _ 

co k = k , g = 


3A ip 2 


cot 


(17) 


(18) 


Note that by working in terms of conformal time and the “conformal field” ip, the frequency 07 , 
equals the (constant) comoving wavenumber k, and does not redshift with increasing a(rf. All of 
the effects of expansion are included in the explicit factors of a(rf and a"(rf. When expansion is 
neglected (i.e., when we set a = 1), equation (|l7|) reduces to the Mathieu equation, the solutions 
of which are exponentially unstable when g < 1 and 7 = 2c+ e, with e < wj,. To see this, 


we may proceed as follows. Our treatment follows the methods outlined in [14], and substantially 
reproduces the results from the alternative approach adopted in Q. 

We introduce the trial solution: 


ipk(v) = c(g)f[a(rf] cos 


Wfc + ) ij 


+ d(rff [a(g)\ sin 


00k + \ e ) r l 


(19) 
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Here we have introduced the ‘scaling’ function /(a) to absorb effects from the expansion of the 
universe. We assume that the coefficients c(rj) and d(rj) are slowly-varying with time, as compared 
with the frequency 7 / 2 , although we make no such assumptions about the behavior of /(a). Then 
it is self-consistent to put d, d! ~ 0(e), and to neglect higher derivatives. In general, resonant 
solutions can be found for frequencies ny /2 ~ u k , but each integer n > 1 corresponds to keeping 
terms of 0(g n ) in the perturbation expansion. [14] Since we are only going to keep terms to 0(g), 
we may make the following approximations: 


cos(77/2) cos(77) = - 003(777/2) + - cos(377/2) ~ - cos (' jq / 2 ), 

sin(77/2) cos(77) = —-sin(77/2) + -sin(377/2) ~ — - sin (777/2). 
Then the equation of motion (|lT|) applied to the trial solution of equation (|l9|) yields: 

u k A(v) sin(7?7/2) + u> k B(rj) cos(77/2) + O [e 2 , g 2 ) = 0, 
with the coefficients A and B given by 


( 20 ) 


( 21 ) 


k 


1 


A — —2c’f — cf ( 2 -|-I -|- (2d! f + df") — df I —gu: k a + e H- 




1 


1 a" 


uj k a 
n 


B — 2d'f + df' ( 2-f-'j -I-(2 c'f' + cf") + cf (-gu k a~ — e -^ . (22) 

V &k) &k V 2 uj k a J 

For equation (21) to be satisfied, we require that both of the coefficients (A, B) of the trigonometric 
terms vanish identically. This gives a pair of coupled differential equations involving the coefficients 
c(r?), d(rj), and f(a). 

For the remainder of this section we will neglect all expansion of the universe, and set f(a) = 


a = 1. Then, if we set c(r?) = Cexp(srj) and d(rj) = Dexp(sr]), equations (|2lD and (|22|) yield 

sC + - gu k + e^j D = 0, 
sD + - (-guk — ej C = 0. 


(23) 


These equations may be solved with 


D 1 

77 = T— = T-t 

C y k 


2 9^k € 

\ \gu k + e ’ 


s = ±- 


1 \n 


gu k - e. . 


(24) 


7 










Thus there exists a growing solution and a decaying solution: 


(v) = CXP J^ [Vk cos(777/2) =F sin( 77 ?/ 2 )] 


nVk 


( 25 ) 


This normalization was chosen, following 0 , so that ~ V’fc Vfc = — 1 - 


From equations (24) and (25) it is clear that the solutions "i/’fc will be exponentially unstable 
whenever s is real. Thus, the parametric resonance will occur only within a small frequency band: 


1 

M < -j9u k - 


(26) 


It will be convenient to introduce the variable £ as 


£ = 


2e 

g^k'' 


(27) 


so that resonance occurs when — 1 < £ < 1. In terms of £, s and yk may be rewritten 

9^k ft - 72 1 ' " 

s = —Vl - £ z , Vk = 


1 - 


(28) 


Note that yk —> 1 as ^ —>• 0, near the center of the resonance band. We further introduce the 
function S±(r 7 ) as 

Si fa) = Vk cos ( 777 / 2 ) =F sin( 777 / 2 ), (29) 


so that near the center of the resonance band, as £ —► 0 , 


|-S'±(77)| 2 -7 1 T sin( 777 ), 


or, averaging over a few oscillations near £ ~ 0 , 


(30) 


<is±(i7)r 


1 . 


(31) 


Near the center of the resonance band, then, the solutions ^( 77 ) become 

(32) 

where the irrelevant phase from the S± term has been dropped. 

We can now make use of a Bogolyuobov transformation to solve for the number of particles per 
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mode, Nl es , produced by the decaying inflaton field during its resonant oscillations. In section 2, we 
made the approximation of an instantaneous phase transition at some time 77 *, which can only be 
appropriate for long-wavelength modes. In the present case, however, the resonant modes cannot 
have arbitrarily long wavelengths: from equations fll8| ) and (^), we require k > 2\e\/g. Instead, 
use may be made of the more elaborate time-dependent Bogolyubov transformation developed in 
Appendix B of |]], the results of which are: 

Nk = \/3 k \ 2 = sinh2 ( sr ?)- ( 33 ) 

As explained in f|, when s is not exactly constant in time, but changes adiabatically (such that 
|s'| <C s 2 ), then this equation for Nk may be modified to 

Nk — sinh 2 ( f s drj \ , (34) 

\J res. band / 

where the integral extends over the resonance band (equation (|26|)). Following [iGj, the (divergent) 
coefficient (1 — i 2 )^ 1 has been dropped, because near \l\ ~ 1, the exponent s no longer changes 
adiabatically. 

Equation ([14|), with s given by equation (p4|), is the main result of this section. Clearly, if 
/ sdr/ is large enough within the resonance band, then the number of particles produced will be 
exponentially greater than the number produced from the non-resonant decay studied in section 
2, given by equation (^). There is an important difference, however, between our result for Nk 
in equation (|||) and the corresponding expression in [6;]. The analysis in [ 16 ) was carried out in 
terms of cosmic time t and the fluctuations 5(f>k] the equations of motion thus resemble equations 
0) and (pj) above, with ojk — k/a, and with the assumption ojk H. As studied in Q, then, the 
resonant modes exit the resonance band due both to the redshift of the physical wavenumbers k/a, 
as well as to the back-reaction from produced particles. This is to be contrasted with the foregoing 
analysis: In our case, by working in terms of conformal time and the “conformal field” ipk, the only 
time-dependence of the exponent s comes from the slow time-dependence of the decaying amplitude 
c p 2 . That is, in our case it is g which is slowly changing with time, and not uJk and it is this change 
of g which causes various modes i[)// to slide outside of the resonance band. The time-dependence 
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of g is further addressed in the Appendix. Meanwhile, we will assume for the remainder of this 
paper that y? 2 (r ]), and hence g{rj) and s(rj), may be treated as slowly-varying with time. With this 
assumption, we may now study what happens to the resonant modes when the expansion of the 
universe is included. 


4 Resonant Inflaton Decay in an Expanding Background 


We return to the full equation of motion for the fluctuations dki equation (0), with u k and g as 
defined in equation (^). We introduce the same trial solution, equation (|I9|), and keep terms to 
0(e, g). This leads to equation (21), with the full coefficients A(jq) and B{g) given by equation f[22|). 
We saw in the previous section that a consistent resonant solution may be found when /(a) and 
a are both set equal to unity. In this section, we keep these terms explicit and general, and show 
that consistent resonant solutions may be found when the expansion of the universe is included. 

We split the coefficients A and B into two terms each: A = A\ + A 2 , and B = B\ + B 2 , with 

Ai = - 2 c f - df , 

= -cf ( 2 + ^ (2d7' + df") - df (a 2 - l) + 


Ao = 


J_af 

ujk a 


B 1 = 2d!f + cf ( -gu k - e 


B 2 = df 


' (2 ^ - \ -\ - {2df' + cf") + cf (-guk (a 2 — l) 

V u k ) Wfc \2 ^ > 


1 a" 

COfo CL 


(35) 


In order for equation (21) to be satisfied, we again require that the coefficients vanish: A = B = 0. 
Because we are interested in resonant behavior, we will further look for solutions when we set each 
of the terms A\, A 2 , B 1 , and B 2 separately equal to zero. Then the coefficients A\ and B\ reduce 
to the same pair of linear differential equations as in equation (E3), which may be solved with the 


same y k and s as given in equation (J24Q . The task now is to use the equations 

(£)♦-(£)]° 


_ CI '4)( 2+ j.) + a 

J J V ^k / ^k . 

J J \ UkJ U k L 
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to find a solution for the expansion factor, f[a(j])]. 

The first step is to notice that several terms in equation ( |36| ) may be combined: 

1 d 




(<V), 


( 37 ) 


(fd) dr] 

and similarly for the corresponding cf -term in the second line of equation (^). Then, because 
the coefficients c(?y) and d(i]) are determined by equation (23), the two equations in ([h]) may be 
combined to give 

d r 


dr] - 


e 2sv f 


a(p) rf 2 + ]-g^l (l - a 2 (??)) e 2sv f(g) 


= E 2 (r,)e 2 s «f(r,). 


(38) 


In the second line, we have used the definitions of y^ in equation (|24| ) and t in equation (27), and 
have introduced the coefficient a(p). For a cosmic-time scale factor a(t) oc t p , or, equivalently, a 
conformal-time scale factor a(g) = {a 0 g) p ^ l ^ p \ the coefficient a(p ) is defined as: 


a" p{2p - 1 ) _2 _ ( \ - 2 

7 = FF' = {p)v ' 

Note that a(p) > 0 for p > 1/2, and that for de Sitter expansion, a(oo) = 2. 

From equation (^), we may write the second-order differential equation for /: 

/" + 2s/' - E 2 (j])f = 0. 


(39) 


(40) 


If we now use the fact that s ~ 0(g, e), and that g <C 1, then we may approximate)] E 2 (g) ~ 
a(p) g~ 2 , and 

f" - a(p) r]~ 2 f ~ 0. (41) 


This equation may be solved with the ansatz f(g) = f3 g Pl ’ 2 , with 

1 


Mi, 2 = 


1 ± ^1 + 4a (p) 


(42) 


1 This approximation for E 2 (r)) is further justified far from the edge of the resonance band, when \l\ 1. Also, 

we may assume, as in section 2, that the phase transition occurs near a(r / t ) = 1. 
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The two possible solutions for p, labeled by subscript 1 and 2, correspond to the choice of ± in 
equation (|42|); the choice of whether pi or p 2 is appropriate will be determined below. With this 
solution for /(?]), the resonant solutions ipk may be written: 


,±, \ f(v) exp[± 577 ] , 

i’k ( ? ?) =- 1= — s ±(v), 


nvk 


(43) 


with S±(rj) given by equation 

We may now choose the appropriate exponent /q _ 2 for the function / by matching the 77 - 
dependence of this resonant solution for t/t/ with the non-resonant solutions found in section 2. For 
fixed 77 , the (non-resonant) long-wavelength solutions behave as (see equations ( 10 ) and (|l 2 |)): 


C oc 


In the resonant case considered in this section, i\) k oc / oc r ] fi1 ’ 2 . 


equation (39), the two exponents may be written 

1 - 2 p 


Ail = 


1 — p 


1 Ai2 — 


p 


1 — p 


(44) 

From the definition of a(p) in 

(45) 


Thus the appropriate exponent for the resonant case is p 2 - Note that _/( 77 ) has the same 77 - 
dependence as the conformal-time scale factor, 0 ( 77 ). Because we require that f[a(r/)] = 1 when 
0 ( 77 ) = 1, we may fix the normalization of / as / = a(rj) = (a 0 'ij) tl2 . In addition, given that a(p) > 0 
for p > 1 / 2 , we will write this exponent as 


M 2 = - 1/421 



1 + 4a (p) 



(46) 


and thus / = (a D 77 )~^ 2 L 

With the ‘scaling’ function / now determined, the full resonant solutions in an expanding 
background may be written: 




exp [±£± 77 ] 



£±(7?), 


(47) 


with the exponent x±(jj) defined by: 


x±(rj) = ±s 


M , / \ 

-ln(a 0 77 ). 

T] 


(48) 


12 











The s-terrn is given by equation (24). As in the previous section, resonance occurs when 2|e| < gtVk- 
With the fluctuations now in the form of equation (|47|), we may again use the result of 

Appendix B of pi], and solve for the number of ^-particles produced per mode k: 


Nk ~ sinh 2 I / x + (rf) dr] 

\J res. band 


(49) 

From the form of g 2 in equation (|46|), it is clear that the expansion of the universe decreases the 
exponential production: f x+drj over the resonance band is less than J sdr), corresponding to the 
non-expanding case of section 3. Nevertheless, it is possible analytically to find self-consistent 
resonant solutions for the fluctuations when the expansion of the universe is included. 

5 Conclusion 

We have demonstrated that the resonance effects studied in [ ( 5] [6] [7] can indeed have dramatic 
consequences for post-inflationary reheating. In particular, it has been shown that solutions for 
the number of particles produced from the resonant decay of the inflaton field may be found given 
a general background expansion, a(t) oc t p . Although the total number of particles produced in 
this expanding case, in equation ([49|), is less than the corresponding expression when H ~ 0 


(equation 34), there still remain regions of parameter space for which this resonant production 
far outweighs the non-resonant production described in section 2. The potential trouble for the 
new resonant reheating scenario is therefore not expansion of the universe, but rather the decaying 
amplitude of the classical 93 -field. This question is discussed further below, in the Appendix. 

As pointed out in || , the case of inflaton decay into further inflatons (as studied above) may be 
of interest for dark matter searches. Following their production, the inflatons would decouple from 
the rest of matter. If the inflatons were given a tiny mass (which has been neglected here), then 
these bosons could serve as a natural candidate for the missing dark matter. Furthermore, much 
of the formalism developed in this paper may be taken over, unchanged, for studying the case of 
inflaton decay into some distinct light species of boson, as studied in (5| [h]. In such decays, it is 
again expected that the parametric resonance would survive a non-zero expansion of the universe. 
The difference for these decays would reside solely in the form of ip 2 (g ). 
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The foregoing analysis for the simple model of a minimally-coupled scalar field cj) with an 
Einstein-Hilbert gravitational action can also be applied directly to many classes of Generalized 
Einstein Theory inflationary models, such as those studied in j22|. For these models, the confor¬ 
mal transformation factor P(x) —> 1 as the inflaton field reaches the minimum of its potential. 
During the ^-field’s oscillatory phase, the non-minimal (j) 2 R coupling takes the form of the stan¬ 
dard (167rG) _1 i? gravitational action of a minimally-coupled theory. Near the epoch of reheating, 
then, the effective lagrangian takes the form of equation H), from which the reheating analysis 
may proceed as above. It would be interesting to calculate Nf~ for the new models of open in¬ 
flation [^], following the methods developed above, which might display qualitatively different 
reheating scenarios. This is the subject of further study. 

6 Appendix 


It is important to understand how ip 2 , and hence g and s, changes with time. A proper treat¬ 
ment would be to extend the Minkowski-spacetime approach adopted in Q |||, and to work out 
the thermal Green’s functions for the decaying inflaton field for the general background expansion 


considered here. (See also [24].) As a first approximation, we may posit the reasonable phenomeno¬ 
logical ansatz, which may be compared with figure 1 in [^]: 

<P 2 {v) - ( P 2 0 e~ r,/T , (50) 

where r is some damping time-scale, most likely to be determined from numerical integration. In 
other words, we will assume here that the frequency of the cos(jrj) oscillations may be treated as 
nearly constant, while the amplitude of the oscillations decreases. From figure 1 in |7|, it is clear 
that we may assume that r~ 1 <C u>k- If we then write g{rj) = g Q exp(—r//r), and keep terms to first 


order in g Q , (e/w*,), and (r _1 /ujk), we may repeat the above analysis in equations (17)-(25). We 
will restrict attention to the non-expanding case, since the expanding case examined in section 4 
reveals the same dependence of Nk on g{rj). 

We introduce the trial solution 


' i /’fc( r ?) = c(rj)e r, / T cos( 7 r// 2 ) + d(rj)e sux^g/2). 


(51) 


14 





Again taking d. d! ~ 0(e) leads to the coupled differential equations: 


- 2d + 2 ct 1 — de — -dg 0 u ) k e V ^ T 
2d' — 2dr~ 1 — ce + -cg o i0 k e~ v ^ T 


0, 


0. 


(52) 


These equations may be solved approximately by writing c{rf) = C exp [0rj\ and d{rf) = D exp [6r)\, 
with 


D 

7=7 = =F- 


1 


C ^y k (v) T \ exp [-77/7] + e’ 

<%) = e - 2? ?/ T -e 2 + i = ±s(r/) + 

With these solutions, the growing and decaying mode functions take the form: 

V’fcfa) = [y fc (7y) cos( 777 / 2 ) =F sin( 777 / 2 )] . 

Vdykiv) 


iff 0 w fc exp[-77/r] - e 


(53) 


(54) 


These solutions have the benefit over the corresponding solutions in equation (25) of ending the 
resonance in a natural way: unlike the solutions from section 3, these mode functions do not 
continue to grow exponentially forever. Now the resonance band is explicitly time-dependent: 


e < 


9o^k c - v /T 

2 


(55) 


In addition, the slight time-dependence of 777,(77) and s{rj) still allows these mode functions to be 
normalized as: V’fcV’fc - d k 'dt = —1 + O (r~ l /u k ). 

With this simple ansatz for the decaying inflaton field amplitude, the exponent s(rj) changes in 
time as 

IVI T~^ 1 

(56) 


s 2 s(rj) [1 — £ 2 (g)]' 

In analogy with the situation in section 3, we have defined the variable £( 77 ) = (. 2e)/(g(r])u k )• 
In order to satisfy the adiabaticity requirement for 5(77) (away from the edges of the resonance 
band), we must require r _1 < g 0 u k - Note that this is more stringent than the original assumption, 
t ” 1 <C uj k . If this new constraint may be satisfied, then the number of ^-particles produced may 
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be calculated from: 


[ s(v) dr l = -y f 

-/res. band ^ J— 1 


re 

y 


+1 d £^H 

£ 2 


— arcsin {£) 


e=+i 


e =-1 


7rre 


(57) 

(58) 


(59) 


This gives N k — sinh 2 (/ sdrf) = sinh 2 (7rre/2). Given the constraints on both r and e, we find 

N k > 1. 

As emphasized in Q [||] , the decay of the amplitude of the inflaton might instead take the form 
of a power-law in time, rather than the exponential form adopted in equation (^). If we write 

^ ^ (it ) ’ 

with b > 0 a real constant, and thus 

g(v)=9o(^)~ b , = 

\VoJ go^k \VoJ 

then we may again repeat the steps in equations (ll7|)-(|25|) . The solutions y k and s again take the 

form 

Vk{v) = 


(60) 


\ 


\g(ri)u k +e 1 If 1 

W vT' s(, ) = T? (, )“‘ 


— e- 


(61) 


and the growing and decaying mode functions again take the form of equation (54), with y k and s 


now given by equation (61). Now the exponent s(rj) changes in time as 

|s'| 6r/ _1 1 


(62) 


s 2 s(rj ) [1 — £ 2 (rj)} 

In order to satisfy the adiabaticity condition, we therefore require b < g Q . Then the number of 
</>-particles produced may be calculated from 

i/& f+i . 


/res. band 


^ - f(y) 7> i/i 


£ 2 


* o(i)-^y 1/b , 


(63) 


where £ 0 = £(rj 0 )- From equation ([K]), we see that £ 0 < ^ max = 1, so that given 0 < b < g Q <C 1, 
£o l ^ b 1. This yields N k ~ sinh 2 (e77 0 £ 0 1 ^/(2 b)) ^$> 1. 
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These two different forms for the decay of the inflaton field due to back-reaction from produced 

, both 


particles, the exponential decay of equation (5C), and the power-law decay of equation 
allow a large production of particles from the inflaton’s decay, so long as the adiabaticity require¬ 
ment on s(rj) may be met. It is clear, however, that a more extended treatment of the effects of 
back-reaction on the decaying amplitude tp 2 {ij) is required, and is the subject of further study. 
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